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We introduce an iterative importance truncation scheme which aims at reducing the dimension of the model 
space of configuration interaction approaches by an a priori selection of the physically most relevant basis states. 
Using an importance measure derived from multiconfigurational perturbation theory in combination with an im- 
portance threshold, we construct a model space optimized for the description of individual eigenstates of a given 
Hamiltonian. We discuss in detail various technical aspects and refinements of the importance truncation, such 
as perturbative corrections for excluded basis states, threshold extrapolation techniques, and different iterative 
model-space update schemes. We apply the idea of the importance truncation in the context of the no-core shell 
model (NCSM) for the ab initio description of nuclear ground states. In a series of benchmark calculations for 
closed- and open-shell nuclei up to '*0 we compare the ground-state energies obtained in the importance trun- 
cated NCSM to the full NCSM. All calculations show an excellent agreement of importance truncated and full 
NCSM for all cases where the latter is feasible. The results demonstrate that the importance truncated NCSM, 
while preserving most of the advantages of the full NCSM, gives access to much larger Ny^JiQ. spaces and 
heavier nuclei. In this way we are able to perform importance truncated NCSM calculations for nuclei like '^C 
and "^O up to A'„,„ = 22. 
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PACS numbers: 21.60.De, 21.60.Cs, 21.30.Fe, 02.70.-c 
I. INTRODUCTION 

Configuration interaction (CI) approaches play an impor- 
tant role for the description of quantum many -body systems in 
many different areas of modern physics, ranging from atomic 
and molecular physics and quantum chemistry to condensed 
matter and nuclear physics. Well known examples for Cl-type 
methods from the different fields include the full and trun- 
cated configuration interaction methods for the many-electron 
problem in molecular physics and quantum chemistry fll], the 
exact diagonalization approaches for Heisenberg- or Hubbard- 
type problems in condensed matter theory Jj-lstl, or the diag- 
onalization shell model or general configuration-mixing ap- 
proaches in nuclear structure physics Hlsll- 

The basic framework of all of these methods is the same: 
Within a model space spanned by a set of many-body states, 
the eigenstates of the Hamiltonian are determined through a 
large-scale numerical solution of the matrix eigenvalue prob- 
lem. The many-body states forming the basis of the model 
space are often Slater determinants of a set of single-particle 
states. The basic parameter which determines the difficulty 
and computational cost of such calculations is the dimension 
D of the many-body model space, i.e., the linear size of the 
Hamilton matrix. If the full eigenspectrum is required, then 
exact numerical diagonalizations are routinely performed for 
dimensions up to D ~ 10^ nowadays. Often only a few eigen- 
states are of interest, such that Lanczos-type algorithms pro- 
vide a very efficient tool and expand the domain of tractable 
model-space dimensions io D ~ 10^ iSEtl ™d possibly 10'° 
through massive parallelization i^. 

We consider applications which require only one or few 
low-lying eigenstates. In those cases the model space often 
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contains a significant number of basis states that contribute to 
the basis expansion of the target eigenstates with extremely 
small or vanishing amplitudes. If these basis states would be 
omitted from the outset, the target eigenstates obtained by a 
solution of the eigenvalue problem in the truncated space and 
all observables derived from them would change only little. 
The diagonalization in the truncated space gives a variational 
approximation to the full eigenstates whose quality is directly 
controlled by the threshold on the amplitudes used to identify 
the important basis states. In order to exploit this idea, we 
need a way to estimate the amplitudes of the individual basis 
states without actually solving the full eigenvalue problem. 
This can be done in the framework of many-body perturba- 
tion theory, using the amplitudes for the first-order perturba- 
tive correction of an initial approximation for the target states 
as an importance measure 13, HD. This is the concept of the 
importance truncation scheme discussed in this paper. 

We will focus on the nuclear many-body problem in the 
framework of a large-scale shell-model approach. However, 
all of the conceptual developments are generic and can be 
applied to any Cl-type many-body method for other quan- 
tum systems as well. For the nuclear many-body problem, 
we aim at an exact ab initio solution for a Hamiltonian in- 
cluding a realistic nuclear interaction. In this context, the 
no-core shell model (NCSM) is the most successful Cl-type 
method at present H HI [H [H [il • The model space 
of the NCSM is spanned by Slater determinants constructed 
from harmonic-oscillator single-particle states with an upper 
limit on the unperturbed excitation energy of the many-body 
basis states of Nmayfi^- A unique advantage of the Nmsufi^ 
truncation is the possibility to separate intrinsic and center- 
of-mass degrees of freedom and thus to obtain translationally 
invariant intrinsic states. The NCSM is able to provide a com- 
plete description of the ground and low-lying excited states in- 
cluding all relevant observables, like energies, transition ma- 
trix elements, form-factors and densities. It has been applied 
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very successfully to nuclei up to mass A ~ 13 using realistic 
Hamiltonians involving two- and three-nucleon interactions, 
e.g. the modern interactions derived within chiral effective 
field theory ill[il[il]. A similar set of observables and nu- 
clei is accessible in Green's function Monte Carlo calculations 
which, however, are restricted to certain classes 
of local interactions. Coupled-cluster methods, which have 
recently been used in connection with chiral two-nucleon in- 
teractions, have provided predictions also for heavier closed- 
shell nuclei ll20l|2lll . 

The range of applicability of the NCSM is limited solely 
by the combinatorial growth of the model space with particle 
number A and energy truncation Nmaxfi^- For '^O the model 
space dimension reaches the order 10^ already for A^max - 8, 
which is typically not sufficient to obtain results that are con- 
verged with respect to A^max- Since the model-space dimension 
is its only crucial limitation, the NCSM provides the optimal 
framework for implementing the importance truncation idea 
As we will discuss in detail, the importance truncated no- 
core shell model (IT-NCSM) obtained in this way extends the 
NCSM to a much larger domain in A and A^max- 

This paper is organized as follows: In Sec. |II]we discuss 
the general elements of the importance truncation scheme that 
can be employed in any Cl-type calculation. In Sec. Hill we 
combine these elements with the NCSM and discuss the basic 
properties of the IT-NCSM. In Sees. llVlandrvjwe present a se- 
ries of large-scale benchmark calculations in the IT-NCSM for 
ground states of different closed- and open-shell nuclei up to 
"'O and compare to the results of the full NCSM. Throughout 
this work we restrict ourselves to a regime where full NCSM 
calculation are still possible to some extent so that a detailed 
assessment of the importance truncation is possible. 



II. IMPORTANCE TRUNCATION SCHEME 
A. Concept 

Consider a quantum many-body system whose ground and 
excited states shall be determined by solving the eigenvalue 
problem of the Hamiltonian in a large model space. The nu- 
clear shell model is a typical example: The model space is 
spanned by a set of Slater determinants of harmonic oscilla- 
tor single-particle states and the lowest few eigenvalues and 
the corresponding eigenvectors of the Hamilton matrix are de- 
termined. Similar configuration interaction (CI) methods are 
used throughout many fields of physics and chemistry. 

In all of these methods the many-body model space is con- 
structed in a combinatorial fashion with some global trunca- 
tion. In the no-core shell model in nuclear physics the model 
space is spanned by all possible Slater determinants con- 
structed from harmonic oscillator single-particle states with 
total excitation energies up to Nmaxft^- In a full configuration 
interaction calculation in quantum chemistry the model space 
is spanned by all Slater determinants that can be constructed 
from a given finite set of single-particle orbitals. 

These global truncations do not account for the specific fea- 
tures of the Hamiltonian and the physical properties of the 



state one is interested in. As a result the model space con- 
tains a substantial number of basis states which are irrelevant 
for the description of a specific eigenstate, e.g., the ground 
state. The basic goal of the importance truncation scheme is 
to identify the important configurations for the description of 
one or a set of target states using the information provided 
by the Hamiltonian. Only the important states are selected 
to construct a new, greatly reduced model-space in which the 
eigenvalue problem is eventually solved. These importance 
selection ideas have been pioneered in quantum chemistry 
in the 1970s leading to a number of diff'erent computational 
schemes (see Sec. Ill Gb . The crucial ingredient is an a pri- 
ori measure for the importance of individual basis states. One 
possible framework to construct a simple importance measure 
is low-order multi-reference or multiconfigurational perturba- 
tion theory as discussed in the following section. Though the 
following is applicable to all types of configuration interaction 
approaches, we will employ the language of the nuclear shell 
model for convenience. 



B. Multiconfigurational Perturbation Theory 

We start from a full model space Alfuii spanned by a set 
of many-body basis states lO,,) — for example the harmonic- 
oscillator Slater determinants of the shell model with some 
model-space truncation. Furthermore we assume a reference 
state l^rei) being a zeroth-order approximation for the eigen- 
state of the Hamiltonian we are interested in, e.g., the ground 
state. In general, the reference state can be a superposition of 
basis states from a subspace A1,ef of the full model space 

veAlrcf 

This initial approximation can be obtained, e.g., from a pre- 
vious CI calculation for a smaller space. In the simplest case 
the reference space can be one-dimensional and the reference 
state is given by a single basis state IOq) corresponding, e.g., 
to the ground state of a closed shell nucleus in an independent- 
particle shell model. 

Now we would like to use many-body perturbation theory 
to estimate the leading corrections to the reference state \^ref} 
resulting from states outside of the reference space. Formally 
this requires the use of multireference or multiconfigurational 
perturbation theory (MCPT) as it is widely applied in quantum 
chemistry 

For setting up the perturbation series we have to split the 
full Hamiltonian H into an unperturbed part Hq and a per- 
turbation W. Since we want to start from the reference state 
I^Pref) as an unperturbed state, the unperturbed Hamiltonian 
has to be chosen such that 

//0l*l'ief) = frefl*l',-ef> (2) 

with an eigenvalue e,ef given by the expectation value with the 
full Hamiltonian H 

er.f^{%-.f\H\^',,i}. (3) 
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Formally, we can write the unperturbed Hamiltonian which 
satisfies the eigenvalue relation as 



(4) 



For simplicity, contributions from states within Aliet which 
are orthogonal to |*Pref ) have been left out, since they will not 
contribute later on. 

The unperturbed energies ey for basis states outside of the 
reference space Alref which enter into the definition of the un- 
perturbed Hamiltonian (|4|i can be chosen freely. This choice 
of the unperturbed energies — and thus of the partitioning of 
the Hamiltonian — has an impact on the convergence behav- 
ior of the perturbation series and a number of different pos- 
sibilities have been studied in this respect ll22ll . In the sim- 
plest M0ller-Plesset-type formulation of MCPT the unper- 
turbed energies are defined as 



ev = fi-ef + Afv 



(5) 



where Afy is the excitation energy of the basis state |Oy) com- 
puted at the level of the independent-particle picture, i.e. us- 
ing the single-particle energies of the underlying basis. When 
using a harmonic-oscillator basis, the single-particle energies 
are just the harmonic-oscillator energies eg = fiQ.{2na + la + 
3/2). When working with a Hartree-Fock single-particle ba- 
sis, these are the Hartree-Fock single-particle energies. 

Alternatively, in an Epstein-Nesbet partitioning, the unper- 
turbed energies of states outside of the reference space are 
defined via the expectation value of the full Hamiltonian 



ey = (Ovl/ZIO,,) , 



(6) 



which appears to be a more natural choice, but does not guar- 
antee better convergence 12211 . For the present application 
computational efficiency is the prime concern, therefore the 
simple M0ller-Plesset-type partitioning ^ is more appropri- 
ate and will be used eventually (cf. Sec. IIII El ). 

Once the unperturbed Hamiltonian is fixed, the perturbation 
W is defined via 
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and we can easily write out the lowest orders of the Rayleigh- 
Schrodinger perturbation series. For the energy the zeroth and 
first-order contributions read 



£("' = (>F,,f|//0|>Pref> = eref 
£*" = <*Fref|H'l^ref> = 



(8) 



as a direct consequence of our definition of the unperturbed 
Hamiltonian. The second-order contribution to the energy as- 
sumes the well known form 
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where we have used that all matrix elements of Hq between 
|*Fref) and the basis states lO,,) i A^ref outside the reference 
space vanish by construction. 

For the many-body states, the zeroth-order contribution is 
just given by the initial reference state 



The first-order correction is given by 



(10) 



<(Dv|//|'y,ef) 



<l>v) 



(11) 



In all of these expressions we can insert the expansion ([T) of 
the reference state I^Pref) in terms of the basis states. Obvi- 
ously, all these relations reduce to ordinary many-body per- 
turbation theory when dealing which a reference state that is 
given by a single basis state, i.e. for l^i-ef) = |<to). 

In the following, MCPT serves two important purposes: (/) 
It provides an efficient way to assess the importance of in- 
dividual basis states outside of the reference space Alief and 
will thus be the main ingredient in the importance truncation 
scheme. (//) It allows for a direct computation of corrections 
to the energy obtained by an initial shell-model calculation in 
a limited reference space Alief, induced by states outside of 
this simple space. 



C. Perturbative Importance Measure 

The central element of the importance truncation scheme 
is an a priori measure for the relevance of individual basis 
states |<l>v) for the description of a specific eigenstate of the 
Hamiltonian. The target state is represented by an initial ap- 
proximation, the reference state I'I'ref), that carries the correct 
quantum numbers. Based on this reference state, multicon- 
figurational perturbation theory provides a natural framework 
for assessing the importance of basis states outside of the ref- 
erence space Alref- 

A simple yet efficient importance measure can be con- 
structed from the expression ( fTTT i for the lowest-order correc- 
tion to the unperturbed, i.e., reference state |*I'ief). The am- 
plitudes of the individual basis states jOy) i Alref in the per- 
turbative correction ( fTTT ) provide a dimensionless measure for 
the relevance of those states. Thus we can use the perturbative 
amplitudes to define an a priori importance measure: 



£v - ^ref 



1-ef) 



mH\%) 



(12) 



Only those basis states with an importance measure |a:v| larger 
than a threshold value k^^^^ are included in the importance- 
truncated model space. This space is tailored for an optimal 
description of the target state for the given Hamiltonian. In 
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contrast to truncation schemes based on global energy cuts, 
the importance truncation criterion is directly governed by the 
Hamiltonian and the target state. The importance threshold 
/fmin controls the size of the model space and will later on be 
varied to investigate the dependence of the observables on the 
truncation. 

By construction, the importance measure Ky characterizes 
the basis states with regard to their relevance for the descrip- 
tion of the eigenstate. This is not the only possible choice. 
One can define a corresponding importance measure for iden- 
tifying the basis states which are most relevant for the descrip- 
tion of the energy. Using the contributions of the individual 
basis states to the lowest-order correction to the energy (|9]l we 
can define the energy-based importance measure 

Since we are aiming at an optimum approximation to the 
eigenstate, which is then used for computing various observ- 
ables other than the energy, the state-based importance mea- 
sure Ky is conceptually superior and will be used in the fol- 
lowing. In practice both measures lead to very similar results 
though the dimensionless state-based importance measure is 
easier to handle (cf. Sec. IIIIEI ). 

It is important to note that for a two-body Hamiltonian the 
importance weight Ky (as well as ^y) vanishes whenever the 
basis state |<t>,,) differs from all of the states in the reference 
space by more than two single-particle states. If we start from 
a single Slater-determinant as reference state, I'I'ref), then only 
\p\h and 2p2/!-excited states with respect to this determi- 
nant can yield non-zero matrix elements for H and thus non- 
vanishing Ky. In order to access 3/?3/i and 4p4/i-excited states 
directly, the second-order perturbative corrections to the am- 
plitude would have to be used. This shows that the construc- 
tion of the importance truncation via perturbation theory nat- 
urally entails a hierarchy of npnh states. Only \p\h and 2p2h 
excitations of I^Fiet) contribute to the leading-order correc- 
tion, 3/73/1 and ApAh excited states first appear in the next- 
to-leading-order, and so on. In order to avoid the computa- 
tionally demanding evaluation of higher-orders of perturba- 
tion theory we embed the first-order importance measure ( fT2t 
into an iterative scheme for the construction of the importance 
truncated space as discussed in Sec. IIIDI 

Although we focus on the description of a single eigenstate, 
the concept of the importance truncation can easily be gener- 
alized to the simultaneous description of several eigenstates. 
Starting from a set of a few reference states l^j"f), we con- 
struct separate importance measures a-*"* for each reference 
state. The corresponding basis state \(i>y) is included into the 
importance truncated space if one of the importance measures 
/fv"' exceeds the threshold /fmin, i.e. if the basis state contributes 
with a sizable amplitude to at least one of the target states. In 
this way, we obtain a model space tailored for the simultane- 
ous description of all target states. 



D. Iterative Model-Space Construction 

Since the importance measure ( fT2] i constructed within 
lowest-order perturbation theory can only be used to extend 
the reference space by \p\h and 2p2h excitations, we adopt 
an iterative procedure to construct the importance truncated 
model space for a given threshold /fmin. Here we discuss a 
simple and universal update scheme applicable for any CI- 
type problem. More specialized update schemes can be de- 
vised for specific models spaces — we will come back to this 
question in the context of the NCSM in Sec. IIII PI 

Assume we start from a single basis state |<l>o) as an ini- 
tial approximation for the target state, e.g., the ground state of 
a closed-shell nucleus. In the first iteration we use this state 
as reference state - IOq) and employ the importance 

measure to construct all \p\h and 2p2h excitations of the ref- 
erence state with \Ky\ > K^^m. Within this new model space 
A1''^(/fmin) consisting of up to 2p2h excitations we solve the 
eigenvalue problem and obtain an improved approximation 
for the target state 

Yj '^v'l'tv) (14) 

veA(['l(/fmi„) 

with amplitudes C{,'^ defined by the eigenvector 

The improved state I^Ft'') obtained in the first iteration is 
used to construct a new reference state for the second 

iteration. In order to accelerate the evaluation of the impor- 
tance measure, we typically do not use the full eigenstate, 
but project onto a reference space Al[.gJ spanned by the basis 
states Yby) e A1'^''(/fmin) with amplitudes Cj,'^ above a refer- 
ence threshold, |Cv''| > Cmin- The new reference state is thus 
defined as 

l*i'I-rf> = <'f Z '^^"l<»v) (15) 

with a normalization constant N'^^^. Typically the reference 
threshold Cmin can be chosen up to 10-times larger than the 
importance threshold /Cmin without affecting the results, we 
will discuss the threshold dependencies in detail later on. As 
in the first iteration, the importance measure is used to con- 
struct all Iplh and 2p2h excitations with |a-v| > k^^^ on top 
of |*I'[.gJ). Since the new reference state already contains up to 
2p2h excitations with respect to the initial Slater determinant, 
the model space AP'^{Kmm) consist of up to ApAh excitations. 
From the solution of the eigenvalue problem we obtain a new 
approximation of the target state 

l^"^>= Yj '^i'''!*^"^ ^^^^ 

yeA1[-l(%,i„) 

with new amplitudes Cy . This improved state again defines a 
new reference state and the previous steps are repeated. 

This scheme is used for a fully adaptive update of the whole 
model space, i.e., in each iteration the importance of all ba- 
sis states is reassessed using the most recent reference state. 
In this way, the impact of the coupling to higher-order npnh 
states is included when selecting states with lower npnh or- 
ders. This relaxation can have sizable effects. 
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E. A posteriori corrections 

Beyond the definition of the importance measure, perturba- 
tion theory can be used to construct a posteriori corrections to 
the CI energies E{K„^i„), which account for contributions from 
basis states that are not included in the importance truncated 
model space AlC/Cmin)- We distinguish two types of correc- 
tions: (;) those accounting for states that were discarded be- 
cause of an importance measure below the threshold and (//) 
those accounting for configurations which would only be gen- 
erated in the next iteration of the update cycle described in 
Sec. IIIDI because of their npnh-OTder. 

An estimate for the energy contribution of basis states with 
non-vanishing importance measure Ia-,,] < K„^i„, i.e. those that 
were excluded from the importance truncated space AiiKmin) 
for given threshold /Cmin, can be obtained from the second- 
order energy correction (|9]l. One can simply add the individual 
energy contributions of the basis states \<^y} i MiK^m)'- 



^excl('''min) 



■ z 

v«A1(/f„,i„) 



KcDvi//i'y.-ef)r 



(17) 



This amounts to adding the energy-importance measures 
defined in Eq. (fT3] l for the excluded configurations. Eval- 
uating this correction during the construction of the impor- 
tance truncated space does not generate any additional com- 
putational effort since the time-consuming matrix element has 
to be computed anyway for the importance measure Ky. 

Generally, the correction AexciC'Cmin) provides only a rough 
estimate for the contribution of excluded states to the energy, 
since only the coupling to the reference state |*Pref) is con- 
sidered but not the coupling to the majority of other basis 
states in MiK^m)- The primary use of this correction relies on 
the formal property that AexciCfmin) has to vanish in the limit 
Kmin " > 0. This makcs it a unique tool for stabilizing the ex- 
trapolation of the CI energy E{Ky^in) to vanishing importance 
threshold — > 0. This constrained threshold extrapolation 
is discussed in detail in Sec. IIII PI 

The effect of higher-order npnh states that would only 
be generated in the next iteration of the model-space up- 
date can also be assessed via the second-order energy cor- 
rection of MCPT. Assume we have performed two iterations 
of the importance-update cycle starting from a single Slater- 
determinant as initial reference state. The importance trun- 
cated space A1'^^'(A-,nin) then contains up to ApAh excitations 
with respect to the initial reference state. In order to estimate 
the effect of 5p5h and 6p6h configurations we can either per- 
form a third iteration to construct Al'^^C/f'min) and solve the CI 
problem or we can apply MCPT on top of the eigenstate l*!*'^^) 
obtained in the second iteration. Based on the second-order 
energy contribution given by Eq. (|9]l, we define the energy 
correction 



Apt = - ^ 



K(D,|//|>y,ef)l' 

- fief 



(18) 



where the reference state is given by the full eigenvector of 
the second iteration, |*I'ref) - I*?'"'), and the sum runs over all 



5p5h and 6p6h configurations. The computational effort for 
evaluating this correction is almost the same as a full iteration 
of the model-space update because of the complexity of the 
reference state. Reference thresholds and extrapolation tech- 
niques can be employed to speed up the calculations also in 
this case. 

Simpler methods for estimating the effects of higher-order 
npnh configurations are used in the context of truncated CI 
calculations in quantum chemistry ^ . Due to their additional 
benefit of restoring size extensivity in truncated CI calcula- 
tions they are commonly referred to as size-extensivity correc- 
tions ll24ll25ll . The simplest class of corrections are the single- 
or multi -reference Davidson corrections ll24ll . which exist in a 
number of different formulations. In the language of quantum 
chemistry, a correction to the energy obtained in the second 
iteration corresponds to a multi-reference situation, where the 
eigenvector of the first iteration |*I''^'') defines the reference 
state and the second iteration includes singles and doubles ex- 
citations on top of this reference state. Out of the different 
forms of multi-reference Davidson (MRD) corrections we use 
the so-called Davidson-Silver or Siegbahn form ll25ll26ll27ll . 
which can be constructed in the context of perturbation theory. 



iMRD 



^E 
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1 -c- 

21 

2C~ - 1 

^'-21 ^ 



(19) 



where £21 - E^^"^ - is the difference of the CI energies 
obtained in the second and the first iteration and 



[2]|2 



(20) 



veMin 



is the total weight with which the configurations in Al''^, i.e. 
those that were already present in the first iteration, contribute 
to the eigenstate after the second iteration. Obviously the eval- 
uation of the MRD correction does not involve any additional 
computational effort. For each value of the importance thresh- 
old A-^in we can extract the correction AMRD('<^niin) using the 
energies and amplitudes of the two last iterations. Eventually 
the MRD correction is also extrapolated to vanishing thresh- 
old /(-„i„ 0. 



F. Properties of the Importance Truncated CI 

Already at this stage we can identify a few general proper- 
ties of the importance truncated CI, which do not depend on 
the details of the physical system or the model space under 
consideration. 

First of all, it is a strictly variational approach. Since we 
determine energies always from a solution of an eigenvalue 
problem of the Hamiltonian in a restricted space, the low- 
est eigenvalue always provides an upper bound for the exact 
ground state energy. Moreover, the Hylleraas-Undheim theo- 
rem ll28ll applies, i.e., the energy of all states is guaranteed to 
drop monotonically with decreasing and is bounded from 
below by the exact eigenvalue E^^^'^^ in the full model space: 



< En(K^m) < E„{k'-) for 



(21) 
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where £'„(/fniin) is the nth energy eigenvalue obtained in the 
importance truncated space Al(/Cn,in). One can view the whole 
importance-truncated CI scheme as a variational calculation 
with an iteratively improved linear trial state. The set of states 
from which the trial state is constructed as a linear superpo- 
sition, is selected using the importance measure based on a 
previous approximation of the target state. 

Second, the iterative construction of the importance trun- 
cated model space will recover the full model space in the 
limit (/fmin, Cmin) after n/2 iterations, where n < A 
is the maximum npnh excitation possible in the full model 
space, when starting with a single basis determinant as initial 
reference. As we will discuss in the context of the NCSM, 
more elaborate choices of the reference state will guarantee 
that this holds even after a single iteration. Together with the 
monotonous behavior of the energy, this limiting property pro- 
vides the foundation for an a posteriori extrapolation of the 
energies for different importance thresholds towards /<•„,]„ — > 0. 

Third, the importance measure (fT2] i is constructed to iden- 
tify states based on their contribution to the expansion of the 
eigenstates and not based on their effect on the energies. Thus 
the importance truncation using Ky is tailored to generate an 
optimal approximation for the eigenstates in a limited model 
space. The energy can be computed from the eigenstates just 
like any other observable of interest. Therefore, from the con- 
ceptual point of view, all observables are accessible with the 
same precision as the energy. 

Finally, an interesting and nontrivial question that was 
raised in Refs. |29, 30] and addressed in detail in Ref. jH] 
concerns the size extensivity of importance-truncated CI cal- 
culations. In simple terms, size extensivity requires that the 
energy obtained in a many -body calculation for a system com- 
posed of two non-interacting subsystems is equal to the sum 
of the energies obtained in separate calculations for the indi- 
vidual subsystems. Whereas full CI is size extensive, a trunca- 
tion of the space at some fixed npnh excitation level destroys 
size extensivity IHIlsl]. Therefore, importance truncated CI 
calculations based on very few iterations of the model-space 
update discussed in Sec. IIIDI can violate size extensivity. As 
discussed in detail in Ref. [8] a computationally simple way 
to restore approximate size extensivity are Davidson-type cor- 
rections as given by Eq. ( fT9] l. In most cases the effect of 
these corrections is small already after two iterative updates 
of the importance truncated space (cf. Sec. IIV Bb . After a 
sufficiently large number of iterations, i.e. once the model- 
space updates have converged, these size-extensivity correc- 
tions ( fT9] l vanish altogether This is in line with the fact, that 
after A/2 iterations at most the importance truncated CI re- 
covers the full model space in the limit (/Cmin, Cmin) ^ and 
thus would be manifestly size extensive. Although the limit of 
vanishing thresholds is realized only through an extrapolation, 
we can nevertheless presume that the importance truncated CI 
provides an approximately size-extensive result after conver- 
gence of the model-space updates and threshold extrapolation, 
simply because it provides an approximation of full CI with- 
out any explicit npnh truncation. 



G. Comparison with other methods 

The idea of an importance selection was pioneered in quan- 
tum chemistry. Already in the late 1960s and early 1970s per- 
turbative importance measures and thresholds were used to 
facilitate large-scale CI calculations ll3ll[32ll. In a set of sem- 
inal papers Buenker and Peyerimhoff |33, 33, HH introduced 
a configuration-selecting multi -reference double-excitation CI 
approach (MRD-CI), which is one of the benchmark methods 
in quantum chemistry up to today. It starts from a multi- 
configurational reference space and adds individual singles 
and doubles excitations employing a selection criterion based 
on the energy-lowering capability of the new configuration. 
The latter can be quantified either by using the perturbative 
second-order energy contribution ( fT3T l or by explicitly eval- 
uating the change of the energy eigenvalue obtained from 
adding the respective configuration. A threshold value on 
this energy-lowering is used to select the important configura- 
tions which are then included in the model space. Already 
in the initial applications of this MRD-CI scheme in Refs. 
ll33[|34i[35ll . powerful threshold extrapolation techniques were 
employed to correct for the effects of excluded configurations 
(cf. Sec, nil Ft . Moreover, size-extensivity corrections as dis- 
cussed in Sec. Ill El can be considered. 

Essentially all conceptual elements of the IT-CI scheme 
are already present in the MRD-CI (although we learned of 
the MRD-CI only after j?! was published). One difference, 
however, lies in the iterative setup we adopt for the IT-CI 
which allows for a systematic improvement of the importance- 
truncated space. Whereas the MRD-CI is typically imple- 
mented as a one-step calculation, the idea of an iterative im- 
provement of the model space has also been used in 
tum chemistry. An example is the CIPSI method 
which uses a CI calculation for a limited model space of 
important configurations and supplements it with a second- 
order perturbative correction for singles and doubles excita- 
tions on top of the CI model space. The CI space is then itera- 
tively enlarged by including those singles and doubles which 
contribute to the first-order perturbed states with amplitudes 
larger than a threshold value. Also this CIPSI scheme con- 
tains many of the relevant ideas employed in the IT-CI. 

Since these early formulations a large number of 
new implementations and variations of the aforementioned 
importance-selection ideas have been developed 111 |39l 
and are being used for the ab initio description of highly cor- 
related problems in quantum chemistry. 

In nuclear physics the use of importance-selection tech- 
niques is not as far developed as in quantum chemistry. How- 
ever, there are some schemes, particularly in the context of 
the valence-space shell model, which employ similar ideas. 
Among those is the Monte-Carlo Shell Model (MCSM) of 
Otsuka et al. ll4lll . It uses the lowering of the energy eigen- 
value caused by adding a test configuration to a set of refer- 
ence states as a criterion for the relevance of this configura- 
tion. However, the crucial element of this method is that the 
test configurations are generated through an imaginary time- 
evolution of the reference set implemented via an auxiliary- 
field Monte Carlo scheme. Due to this stochastic sampling 
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the individual configurations are no simple shell-model basis 
states anymore, but more complex states containing informa- 
tion on the Hamiltonian already. For the final diagonalization, 
typically supplemented by an angular momentum projection, 
a small number of those MCSM configurations is sufficient to 
capture the relevant physics. 

Another importance sampling scheme has been proposed 
by Andreozzi et al. ll42ll in connection with an iterative 
method for the solution of the eigenvalue problem ll43ll . Here 
the approximations of the eigenvalues obtained during the it- 
erative solution are used to apply an energy threshold crite- 
rion to discard irrelevant states. Horoi et al. have devised a 
truncation scheme based on the diagonal matrix elements of 
the Hamiltonian and applied it in sd and fp-shell calculations 
H. 



III. IMPORTANCE TRUNCATED NO-CORE SHELL 
MODEL 

As the primary apphcation we study the importance trun- 
cation scheme in connection with the no-core shell model 
(NCSM) 101. Applications of the importance truncation in nu- 
clear CI approaches based on a different definition of the full 
model space have been presented in Ref. i^. 



A. Model space 

The NCSM is based on an expansion of the many-nucleon 
state in a basis of Slater determinants of harmonic oscilla- 
tor single-particle states. The model space of the full NCSM 
is restricted solely with regard to the maximum number of 
harmonic-oscillator excitation quanta, A^max, in the many-body 
basis state. In other words, all harmonic-oscillator Slater 
determinants with unperturbed excitation energies of up to 
A^max^fi are included in the model space. 

The combination of harmonic oscillator basis and Nmaxfi^ 
truncation has a unique advantage. Only this model space al- 
lows for an exact separation of the center-of-mass and intrin- 
sic component of the many -body state for all A^max- There- 
fore, one can guarantee that the intrinsic part of the state is 
free of spurious center-of-mass contaminations. Any other 
single-particle basis, e.g. a Hartree-Fock basis, or a differ- 
ent model-space truncation, e.g. a truncation at the level of 
the single-particle states like in other CI methods, will destroy 
this property and induce center-of-mass contaminations of the 
eigenstates which can severely affect intrinsic observables. 

The dimension of the Nmaxfi^ model space grows factori- 
ally with A^max and particle number A. Therefore, full NCSM 
calculations are computationally feasible only for relatively 
light nuclei or in very small spaces. Model space dimensions 
of the order of 10^ are used routinely with present NCSM 
codes i^Ult]. For '^O this allows for calculations in an 
space, which for most realistic Hamiltonians is not sufficient 
to reach convergence. The dimension of the IQTiQ. model 
space is larger than 10"^ and thus just beyond the reach of 



the full NCSM at present. For heavier nuclei the situation be- 
comes progressively worse. 

The importance truncation can be used to efficiently reduce 
the dimension of the Nmaxft^ model space to a tractable size. 
Note that the Nm^^hD. space already reflects a simplistic im- 
portance selection of the individual many-body basis states. 
Based on the perturbative arguments of Sec. Ill CI the ampli- 
tudes of basis states with large unperturbed excitation ener- 
gies will be suppressed by the energy denominator in ( fT2b . 
Precisely those states are discarded through the A^niax^^^ trun- 
cation. However, the numerator of ( fTSl i and thus the full 
Hamiltonian, is not considered in this simplified picture. The 
Nnmxfi^ model space is not adapted to the specific properties 
of the Hamiltonian or the target states under consideration. 
By using the importance truncation in combination with the 
NmsLxfi^ model space we also include these aspects. 



B. Hamiltonian 

For the following discussion we use a translationally invari- 
ant Hamiltonian composed of intrinsic kinetic energy Tint - 
T - Tcni and a realistic two-nucleon interaction Vnn: 

21'* ^ 

//int = Tim + ^NN = ^ ^ J] (IJj + ' (22) 

where q,j = i(p,- - p^) is the relative two-body momentum 
operator and fi = mN/2 the reduced mass. 

In principle any two-body interaction can be used as input. 
In this work we restrict ourselves to unitarily transformed in- 
teractions derived in the framework of the Unitary Correlation 
Operator Method (UCOM). Starting from the Argonne V18 
potential a unitary transformation is used to account for short- 
range central and tensor correlations leading to a phase-shift 
equivalent effective interaction with improved convergence 
properties. The conceptual details of the UCOM approach 
are discussed in Refs. [46, 47, 48J. Further details regarding 
the calculation of matrix elements of the Vucom interaction 
and the determinantion of the optimal correlation functions 
are discussed in Ref. ll49ll . 

For all of the following calculations we use the 'standard' 
set of correlation functions introduced in Ref. ll49ll with 
a triplet-even tensor correlator with range parameter I^t = 
0.09 fml This value was chosen such that experimental bind- 
ing energies for ^H and "*He are roughly reproduced in full 
NCSM calculations. Though improved correlation functions 
are available ifsoll . there exists a number of different many- 
body calculations for this first-generation Vucom interaction. 
Calculations for light nuclei in the NCSM and other methods 
ifsill show that Vucom exhibits good convergence properties 
and provides a realistic description of a number of observ- 
ables. Studies of heavier nuclei in Hartree-Fock plus second- 
order many -body perturbation theory demonstrate that this in- 
teraction provides reasonable binding energies throughout the 
whole nuclear mass range without the explicit inclusion of a 
three-body interaction |52]. Therefore the Vucom interaction 
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provides a realistic testbed for the many-body methods inves- 
tigated here. 

We emphasize that all of the following calculations use the 
Hamiltonian ( l22l i without further transformations, i.e., there 
is no additional Lee-Suzuki similarity transformation as in the 
ab initio NCSM S Eil IH H S El . Here the term 
NCSM solely refers to a Cl-type calculation specifically using 
an Nms.}fiQ. model space. 



C. Implementation 

The implementation of an importance truncated NCSM dif- 
fers from a conventional NCSM code. The computationally 
most demanding part is the construction of the importance 
truncated space itself. Due to the reduction of the dimen- 
sion of the model space, the subsequent computation of the 
Hamilton matrix and the solution of the eigenvalue problem 
are simpler than in a full NCSM approach. 

For generating the importance truncated space for a given 
reference state I^Fref), we use an algorithm motivated by the 
structure of the importance measure (fT2l i. We loop over all 
basis states lO^) e A1,-ef contained in the reference state |*P,-ef) 
and create all \p\h and 2p2h excitations of each of them. In 
order to avoid creating duplicates, we discard any newly cre- 
ated determinant lO^) that has a non-vanishing matrix element 
of the Hamiltonian with any of the states in Alief that were 
processed previously. This update scheme, which is also used 
in importance selecting CI approaches in quantum chemistry 
ll33i [36l I40I1 . is much more efficient than the simple scheme 
employed in There, explicit loops over all possible npnh 
excitations of the independent-particle shell-model state |<l>o) 
were used to generate candidate states for evaluating the im- 
portance measure. Though duplicates are excluded from the 
outset, this scheme becomes less efficient already at the 3p3/i 
order and it eventually limited the calculations in Ref. |l7i| to 
states up to the ApAh level. Therefore, all results presented 
here are based on the refined implementation without any ex- 
plicit limitation of the npnh level of the states considered. 

Evidently, the cost for the model-space update grows 
quadratically with the number of basis states in the reference 
I^Pref). Therefore, as discussed in Sec. IIIDI we introduce an 
additional reference threshold and define the reference state 
I^Fief) using the dominant components of the previous eigen- 
state. Typical reference thresholds Cmin are of the order of 
lO""* which leads to reference states composed of typically 10^ 
basis determinants. We always check that a further lowering 
of the reference threshold does not produce sizable effects. 

Eventually we obtain a list of basis states spanning the im- 
portance truncated model space including their importance 
weights. The typical dimensions we deal with are of the 
order 10^. These problems can be handled by conventional 
Lanczos- or Arnoldi-type algorithms — in addition to sim- 
ple Lanczos-implementations we use the implicitly restarted 
Arnoldi algorithm of the ARPACK library pSll. The many- 
body matrix elements of the Hamiltonian are pre-computed 
and stored in memory or on disk. Using the known impor- 
tance weights as initial pivots, one can obtain convergence 



of a single target states after typically 10 iterations. Eventu- 
ally, we obtain energy eigenvalues and amplitudes of the target 
states. Since the eigenstates are — at no additional cost — given 
in a simple shell-model representation, we can easily use them 
for subsequent computation of various expectation values and 
density distributions or form-factors. 

The time-consuming parts to the code, i.e. the construction 
of the importance-truncated space and the computation of the 
Hamilton matrix, can be easily parallelized with practically 
no communication overhead and perfect scaling. We use a 
hybrid OpenMP plus MPI parallelization strategy to make op- 
timal use of the memory resources of modern multi-core ar- 
chitectures. As compare to a typical full NCSM, the particle 
numbers and model space sizes in the importance-truncated 
NCSM are not limited by the available memory. Larger model 
spaces or particle numbers only require more CPU-time for 
the construction of the importance-truncated model space. 



D. Iterative Construction of Model-Space: IT-NCSM(0 vs. 
IT-NCSM(seq) 

We can use the universal update scheme described in 
Sec. IIIDI for the iterative construction of the importance trun- 
cated NmsLxfi^ space for any given A^max- For targeting the 
ground state, we would start with a OHO. eigenstate as the 
initial reference state — for a closed-shell nucleus this is just 
the independent-particle shell-model determinant. In a first 
iteration the importance update is used to generate all rele- 
vant I pill and 2p2h excitations within the Nmaxfi^ space un- 
der consideration. Using the eigenstate in this importance- 
truncated space as reference state, a second iteration will give 
access to all basis states up to the 4p4h level with respect to 
the initial OftQ state. Typically two or three iterations of the 
importance update cycle are sufficient to obtain convergence, 
i.e. a result which is not changed anymore by another impor- 
tance update. In the following we will identify those calcula- 
tions with the label IT-NCSM(;), where / indicates the number 
of iterations. 

However, for the N^nsofi^ space of the NCSM there exists 
a more efficient alternative. Typically we are interested in a 
sequence of calculations for growing A^max in order to assess 
the convergence behavior with increasing model-space size. 
We can combine this sequential increase of A^max with the im- 
portance update in an elegant way. Assume we start with a 
complete NCSM calculation in a OfiO or 2/iO space. Using 
the eigenstate obtained in this small, say 2fi^ space as refer- 
ence state we construct the importance truncated 4/if2 space 
and solve the eigenvalue problem again. The resulting eigen- 
state then defines the reference state for the construction of 
the importance truncated 6MI space, and so on. We will iden- 
tify calculations based on this sequential update scheme by 
IT-NCSM(seq) in the following. 

This sequential scheme has an important conceptual ad- 
vantage; The maximum npnh excitation with respect to the 
O/jQ space that is contained in an Nmaxfi^ space is of or- 
der n = A^max- Therefore in each step of sequence A^max = 
0,2,4,6,... the maximum npnh-oidei increases by 2 and a 
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FIG. 1: (color online) Correlation between the importance measure 
Ky and the amplitude C,, obtained by solving the eigenvalue problem 
in an IT-NCSM(2) calculation of '^O with N^,,^ = 8 and nil = 22 
MeV. The panels correspond to the dilFerent «p«/j-orders as indi- 
cated. 



single importance update at each step is sufficient to access 
all npnh orders that can appear Thus, the sequential up- 
date scheme recovers the complete Nm^xfiQ. model space in 
the limit (/Cmin, Cmin) — * and does not impose any explicit 
limitation regarding the «/9«/i-content of the space. We need 
to apply the importance update only once for each value of 
A^max, in the iterative scheme we would need / = A^max/2 it- 
erations to formally achieve this. We will apply and compare 
both schemes in Sec. HV] 



E. Importance Measure 

As a first test of the reliability of the importance measure, 
we can compare the perturbative estimate Ky for the amplitude 
of a given basis state \Oy} with the amplitude Cy resulting 
from the diagonalization. Whereas the a priori importance 
measure Ky only includes the coupling to the states from the 
reference space, the a posteriori amplitudes Cy are affected by 
the mutual coupling of all states. Nonetheless, the Ky provides 
a reasonable estimate for the amplitudes Cy which is sufficient 
to identify the important basis states. 

This is illustrated in Fig. [1] for an importance-truncated 
NCSM calculation for in an 8^Q space using two itera- 
tions of the importance-update of the model space for an im- 
portance threshold K^i„ - 5 x 10"^. The coiTelation plots re- 
late the importance measure Ky of the individual basis states 
with the corresponding amplitudes C,, in the final eigenstates. 
There is a clear correlation between the two quantities which 
is sufficient to predict which basis states are important for an 
adequate representation of the final eigenstate. The scatter- 




FIG. 2: (color online) Ground-state energy of '*0 obtained in IT- 
NCSM(2) calculations for Mnax = 8 and hO. = 22 MeV as function 
of the dimension D of the importance-truncated space. The different 
symbols correspond to different definitions of the importance mea- 
sure: the state-based importance measure Ky (•), the energy-based 
importance measure Xv (♦). and the state-based importance measure 
defined in an Epstein-Nesbet partitioning (a). 



ing around the diagonal reflects all couplings that are not ac- 
counted for in the lowest-order perturbative estimate. 

As mentioned in Sec. Ill CI there are other options to define 
an importance measure in the framework of multiconfigura- 
tional perturbation theory. A natural alternative to the state- 
based importance measure Ky is the energy-based importance 
measure Xv defined in Eq. ( fT3l l. One could also consider an 
Epstein-Nesbet partitioning as discussed in Sec. lIIBI to set up 
the perturbative corrections and define a state-based impor- 
tance measure k^^. 

In order to assess the efficiency of the three measures we 
perform a series of calculations with different values of the re- 
spective importance thresholds and plot the energy eigenvalue 
versus the dimension of the importance truncated space as a 
parametric curve spanned by the importance thresholds /Cmin, 
Xmin^ and K^^, respectively. Since the whole approach is vari- 
ational, the measure which leads to the lowest ground-state 
energy for a given dimension D of the importance-truncated 
space is most efficient in selecting the D most important basis 
states. 

An example of this analysis is shown in Fig. |2] again for the 
ground state of '^O in an SSQ space. In all cases the NCSM 
ground state in a complete 2HQ. space was used as reference 
state for the construction of the importance-truncated space. 
The points obtained with all three definitions of the impor- 
tance measure essentially fall onto the same line, i.e. all mea- 
sures are able to identify the most important configurations 
with the same efficiency. We therefore use the conceptually 
and computationally simplest importance measure, the state- 
based measure Ky of Eq. ( fT2] i in all following investigations. 



F. Threshold Dependence & Extrapolation 

The variation of the threshold Kn^i„ is an important probe 
for the quality of the importance truncation and the basis for 
an extrapolation to vanishing threshold /Cmin — > as it will be 
used later on. To this end, all IT-NCSM calculations are per- 
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FIG. 3: (color online) Threshold dependence of the energies and the 
model space dimension for a IT-NCSM(2) calculation of '^0 with 
A'max = 8 and TrCl = 22 MeV. (a) Energy eigenvalues and as func- 
tion of without (•) and with (♦) perturbative correction for the 
excluded configurations, (b) Total dimension of the importance trun- 
cated space (•) as well as the number of 2p2h (o), 3p3h {<)), and 
4p4/7-configurations (□) with varying A-^i,,. 



formed for a sequence of different values for K^i^. For each 
threshold value the importance truncated space is different and 
the eigenvalue problem has to be solved again. However, this 
can be done at small computational cost. The importance trun- 
cated space and the Hamilton matrix are initially determined 
for the smallest After the solution of the eigenvalue prob- 
lem for this threshold, all basis states that are not part of space 
for the next-larger importance threshold and the correspond- 
ing matrix elements are removed, and the eigenvalue prob- 
lem is solved again. Hence, the time consuming construction 
of the importance truncated space and the computation of the 
Hamilton matrix is done only once for a whole threshold se- 
quence. 

The dependence of the energy and of the model-space di- 
mension on the importance threshold K^in in IT-NCSM(2) cal- 
culations for ""O with different A^max is illustrated in Fig. [3] 
The energy eigenvalue ^(/Cniin) obtained in the importance 
truncated space decreases monotonically with decreasing K,^in 
as expected from the variational principle and the Hylleraas- 
Undheim theorem. At the same time, the dimension of the 
importance truncated space increases exponentially with de- 
creasing Kmin- The number of configurations of higher npnli- 
order in particular grows rapidly as the threshold is lowered. 
This behavior reflects the mechanism behind the importance 
truncation scheme; The configurations which are most impor- 
tant for the description of the target state have large Ky and are 
included already for large thresholds. With decreasing thresh- 
old ATmin basis states of lesser importance are successively in- 
cluded. Their number increases dramatically but the effect 
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FIG. 4: (color online) Threshold extrapolation of the ground-state 
energy of '^O (JiQ = 22 MeV) obtained in IT-NCSM(2) for differ- 
ent A'max- Shown are the perturbatively corrected energies £,i('<:miii) 
as function of /Cmin for A = 0, 0.5, 1, 1.5, and 2 (data sets from top to 
bottom within each panel). For A = (•) the original energy eigen- 
value ii(A'min) is recovered, for A = 1 (♦) we obtain the perturbatively 
corrected energy £(/<-n,i„) -I- Aexd(fmin)- The lines show the results of a 
simultaneous constrained fit for all data sets using 4th order polyno- 
mials (see text). 



on the state and the energy remains moderate, facilitating ap- 
proximations to estimate their effect on the energy without 
including them explicitly in the model space. 

The simplest approximate way to account for the excluded 
basis states is the a posteriori energy correction Aexci('<'min) 
given by Eq. ( [TTl i on the basis of the second-order MCPT 
contribution. The corrected energies ^(/Cmin) + AexciCfmin) are 
also depicted in Fig. 13 a). Although AexcKA^min) provides only 
a rough estimate for the contribution of excluded states, the 
A-,nin-dependence of the corrected energy is much weaker than 
the dependence of the uncorrected eigenvalues E{Kmm) — if the 
correction were exact we would expect the corrected energies 
to be independent of A-„,in. In many cases the corrected energy 
at a single value of A-„,in can already serve as a good approxi- 
mation for the full results in the limit K^i„ — > 0. 

A more reliable way to recover the contribution of excluded 
configurations is an a posteriori extrapolation of the energies 
to vanishing importance threshold. Due to the smooth and 
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monotonic behavior of the energies E(Kmm) one can attempt a 
direct numerical extrapolation k^[„ — » as done in Ref. |01- 
Since the general shape of the E(Kmi„) curve varies, we will 
generally use polynomials in /Cmin fitted to a sufficiently large 
number of different threshold values for the extrapolation. In- 
stead of ^(A-tnin) one can extrapolate the perturbatively cor- 
rected energy, E(K^in) + Aexci('<'min), which shows a weaker 
threshold dependence than the eigenvalues and, therefore, al- 
lows for a more stable extrapolation. The extrapolation can be 
stabilized further by performing a simultaneous fit of E(Kj„i„) 
and ^(/Cnjin) + AexcK^min)- Sincc the perturbative correction 
Aexci(fmin) has to vanish in the limit /Cmin both extrapo- 
lations should formally give the same value at K^nin = 0, in- 
dependent of the absolute quality of the perturbative estimate. 
The formal property £(0) - E{0) + Aexci(O) is used as a con- 
straint in the simultaneous fit and reduces the uncertainties of 
the threshold extrapolation significantly. 

One can even go one step further and define a family of 
energy curves EAiK^in) = E(K^i„) + /lAexci('fmin) with a con- 
trol parameter A. Independent of the choice of A the formal 
property £^(0) - E{0) holds. Using this as a constraint in 
a simultaneous ;i'^-fit of a set of curves for several values of 
A provides very robust extrapolation results. This technique 
has been pioneered by Buenker and Peyerimhoff in the early 
applications of configuration-selecting CI approaches in quan- 
tum chemistry Isill . It solely relies on the fact that the correc- 
tion AexcK'fmin) IS a monotonous function which goes to zero 
(smoothly) as /Cmin 0. 

Examples for this type of threshold extrapolation in the 
case of 1T-NCSM(2) calculations for '^O in diff^erent N^^JiD. 
model spaces are presented in Fig. HI The starting point are 
the energies £(/Cmin) and perturbative corrections AexcK^min) 
obtained for a sequence of importance thresholds in the range 
from /fmin - 3 X 10"^ to 14 X 10"^. Using this input we 
construct data sets for the corrected energies E^iK^in) for 
A = 0, 0.5, 1, 1.5, and 2 and simultaneously fit each of the sets 
by a 4th order polynomial under the constraint that all curves 
meet at /Cmin = 0. The individual data sets and the polyno- 
mial fits are shown in Fig.|4] It is evident that this extrapola- 
tion scheme is most stable if the curves approach the common 
£^(0) value more or less symmetrically. This is the reason for 
the particular set of /l-values adopted here. 

We employ the following threshold extrapolation protocol 
for the applications presented in Sec.|IVl Using a sequence of 
12 equidistant threshold values in the range /Cmin = 3 x 10"^ to 
14x10"^ we perform a constrained simultaneous fit of the cor- 
rected energies E^iKmin) for a sequence of at least 5 different A- 
values using low-order polynomials. The set of /l-parameters 
is chosen such that the common point of all fit curves at 
Kmin = 0, which givcs the final threshold-extrapolated energy, 
is approached symmetrically. In order to assess the uncer- 
tainty of the extrapolation, we drop the smallest and the largest 
value, respectively, from the /l-sequence and perform the si- 
multaneous fit for the remaining data sets. The variance of 
this set of extrapolations defines an uncertainty interval for 
the threshold extrapolated energy. 

Exceptions are very light nuclei, e.g. "^He, where a direct 
extrapolation of the energy eigenvalue £(/<•,„;„) without using 



the perturbative correction AexdC'^min) provides a more stable 
result. The reason is the /Cmin-dependence of AexciCfmin), which 
in very small spaces shows structures that interfere with the 
polynomial extrapolation. 

IV. APPLICATIONS & BENCHMARKS: MAGIC NUCLEI 

We employ the IT-NCSM now for the series of calcula- 
tions for the 0^ ground state energies of various closed and 
open shell nuclei in the p-shell. The aim is to compare the 
results to the full NCSM in diff'erent cases in order to demon- 
strate the robustness of the importance truncation scheme. All 
full NCSM calculations presented in the following were per- 
formed with the Antoine code ll45ll . 



A. Helium-4 

As the simplest benchmark we study the ground-state en- 
ergy of '^He using the Vucom interaction. In this case full 
NCSM calculations can be performed up to very large Nmaxfi^ 
spaces such that convergence is observed. Furthermore, other 
few-body methods, e.g. the hyperspherical harmonics basis 
expansion ifsill . have been employed and yield an independent 
reference value for the ground-state energy. 

First we consider the simple iterative scheme IT-NCSM(/) 
for the construction of the importance truncated model space. 
For fixed A^max we perform up to three iterations of the im- 
portance update starting with the Slater determinant of the 
independent-particle model as initial reference. In each itera- 
tion we solve the eigenvalue problem for a sequence of impor- 
tance thresholds in the range /Cmin = 3 x 10"^ to 14 x 10"^ and 
extrapolate the eigenvalues EiK^i^) to the limit of vanishing 
threshold K^in — > as discussed in Sec. IIII Fl For very light 
nuclei like '^He the direct extrapolation of ^(/Cmin) without per- 
turbative corrections for excluded configuration provides the 
most stable results. For the definition of the reference state 
for the next iteration a reference threshold Cmin = 5 x 10 is 
used. 

The threshold-extrapolated ground-state energies and the 
dimensions of the maximum importance-truncated model 
spaces as function of N^-dx are depicted in Fig. |5] The con- 
vergence with respect to the importance updates of the model 
space is very fast. After two iterations, i.e. at the IT-NCSM(2) 
level, we already obtain stable results which is within 100 keV 
of the full NCSM result. The third iteration only lowers the 
ground-state energy a little further bringing it into excellent 
agreement with the full NCSM, as seen in Fig. |5jb). In the 
case of '^He this convergence pattern may be expected. Af- 
ter two iterations the importance truncated space contains up 
to 4p4h excitations, i.e., the full model space can be gener- 
ated in the limit of vanishing thresholds. The minimal change 
in the third iteration is due to a relaxation of the importance 
truncated space, i.e., through the reassessment of the impor- 
tance of all basis states with respect to a new reference state, 
which includes all possible npnh-oideis, the importance trun- 
cated space is better adapted . Further importance updates do 
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FIG. 5: (color online) Ground-state energy and model-space dimen- 
sion as function of A'max for "^He obtained within the IT-NCSM(;) 
scheme for ; = 1 (•), / = 2 (♦), and i = 3 iterations (■) using the 
VucoM interaction for fin = 40 MeV. Panels (a) and (b) show the 
ground-state energies on different scales, including the uncertainty 
estimates for the threshold extrapolation. Panel (c) depicts the maxi- 
mum dimension of the importance-truncated model space. For com- 
parison the results of full NCSM calculations for the same Hamilto- 
nian are included (+). 



not change the resulting energies anymore. 

The agreement with the full NCSM demonstrates the ef- 
ficiency of the importance measure and the reliability of the 
threshold extrapolation. The dimension Dmax of the largest 
model space considered for the threshold extrapolation is up 
to two orders of magnitude smaller than the dimension of the 
full NCSM space, as illustrated in Fig. He). Note that the 
full NCSM dimension is obtained by exploiting all relevant 
symmetries, including parity and time-reversal, to reduce the 
dimension of the eigenvalue problem — it corresponds to the 
'effective dimension' used by the Antoine code. Thus this 
substantial reduction of the model-space dimension by the im- 
portance truncation goes beyond generic symmetries and re- 
ally exploits the specific properties of the Hamiltonian. 

As an alternative to the simple iterative model-space up- 
date at fixed A^niax we can perform these calculations using 
the sequential model-space update IT-NCSM(seq) proposed 



FIG. 6: (color online) Ground-state energy, (a) and (b), and model- 
space dimension (c) as function of A'max for ''He obtained within the 
IT-NCSM(seq) scheme (•) using the Vucom interaction for Tifl = 40 
MeV. For comparison the results of full NCSM calculations with the 
same Hamiltonian are included (+). 



in Sec. IIII Dl Starting from the QhQ space we use the im- 
portance measure to construct an importance truncated 2ftD. 
space. This is used as reference space to construct the impor- 
tance truncated 4HQ. space, and so on. As before we use a 
reference threshold of Cmin = 5 x 10 and a sequence of im- 
portance thresholds starting from /Cn,in = 3 x 10"^. The results 
for the ground-state energies of ''He are summarized in Fig.|6] 
and compared to the full NCSM. The IT-NCSM(seq) scheme 
leads to the same excellent agreement with the full NCSM as 
the IT-NCSM(3). However, the IT-NCSM(seq) is computa- 
tionally more efficient, since only one importance update is 
needed for each value of A'max. 

The dependence of the ground-state energy obtained in the 
IT-NCSM(seq) on the oscillator parameter ftQ is depicted in 
Fig. I?] The comparison with the full NCSM results shows that 
the excellent agreement persists for all frequencies hQ.. The 
particular oscillator frequency HQ. = 40 MeV used in Figs. |5] 
and|6]corresponds to the minimum for the larger space. 

In order to compare our results with other many -body meth- 
ods and with experiment, we perform an exponential extrapo- 
lation of the IT-NCSM(seq) energies for JiQ. - 40 MeV. Since 
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FIG. 7: (color online) Ground-state energies of "^He obtained for the 
VucoM interaction as function of the oscillator frequency TiQ. for dif- 
ferent A'lnax^n model spaces. Shown are the results of IT-NCSM(seq) 
calculations (full symbols) in comparison to full NCSM calculations 
(crosses). 



the calculations are practically converged with respect to A^max 
the main purpose of the extrapolation is to smooth out the 
fluctuations due to the uncertainties of the threshold extrap- 
olation. Using the five data points from A^niax = 16 to 24 we 
obtain a "'He ground-state energy of -28.52(10) MeV. This is 
in excellent agreement with the value of -28.57 MeV that was 
obtained previously in the framework of the hyperspherical 
harmonics approach using the same Vucom interaction ifsill . 
The comparison to the experimental binding energy of -28.29 
MeV only reveals the rough nature of the adjustment of the 
UCOM tensor correlator range that was used in Ref. ll49ll 
to fix the Vucom interaction. In principle one could select 
such that the experimental *He binding energy is reproduced 
exactly. 



B. Oxygen-16 

The ground state of '^O poses a more challenging problem. 
At present, full NCSM calculations can be done routinely for 
spaces up to A^max = 8 with an effective dimension of almost 
0.6 X 10^. For A^max = 10 and 12 the effective dimension grows 
to 1.4 X 10'° and 2.4 x 10", respectively, which is clearly be- 
yond the reach of present NCSM codes. The importance trun- 
cation is crucial in this domain and enables us to treat model 
spaces of up to A^max - 22 and beyond. This limit is set by the 
available two-body matrix elements and not by the IT-NCSM 
calculation itself. 

As for ''He, we first consider the simple iterative IT- 
NCSM(0 scheme using up to two iterations for each A^max 
to construct the importance truncated model space. We use 
a set of 12 equidistant importance thresholds in the range 
^■min = 3 X 10"^ to 14 X 10"^ as input for the simultaneous 
threshold extrapolation as discussed in Sec. IIII Fl The refer- 
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FIG. 8: (color online) Ground-state energy and model-space dimen- 
sion as function of A'max for '^O obtained within the IT-NCSM(;) 
scheme for i = 1 (•) and ; = 2 (♦) iterations using the Vucom in- 
teraction for no. = 22 MeV. In addition the IT-NCSM(2) energies 
after inclusion of the multi-reference Davidson correction are shown 
(■). Panels (a) and (b) show the ground-state energies on different 
scales, including the uncertainty estimates for the threshold extrapo- 
lation. Panel (c) depicts the maximum dimension of the importance- 
truncated model space. For comparison the results of full NCSM 
calculations for the same Hamiltonian are included (+). 



ence threshold is set to Cmin =5x10 . 

A summary of the IT-NCSM(0 results for the ground-state 
energies of "'O up to A^max = 18 is presented in Fig. [8] se- 
lected numerical values are given in Tab. |T] As for the much 
lighter nucleus ''He the convergence of the iterative impor- 
tance updates is excellent. Already after two iterations, i.e. 
for IT-NCSM(2), the full NCSM energies up to A^^ax = 8 are 
produced to an absolute accuracy of better than 600 ke V. 

Instead of performing a third iteration expUcitly, we can use 
computationally simpler estimates for the small correction re- 
sulting from 5p5h and 6/76/1 configuration that are not present 
in the IT-NCSM(2) model space. As discussed in Sec. Ill El the 
simplest a posteriori correction is the multi-reference David- 
son correction (MRD) given by Eq. (fT9T l since it does not 
require any additional computation beyond IT-NCSM(2). The 
MRD corrected IT-NCSM(2) energies are also shown in Figs. 
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TABLE I: Ground-state energies (in units of MeV) for O obtained 
for the VucoM interaction at TiSl = 22MeV with different levels of 
the IT-NCSM. For the IT-NCSM(0 results for i = 1 and 2 itera- 
tions are shown. Furthermore IT-NCSM(2) results with the MRD 
correction l |19t and the perturbative correction J18t for the effect of 
the next iteration are reported. For the IT-NCSM(seq) two different 
reference thresholds have been used: (a) Cmin = 5 x 10""* and (b) 
Cmin = 3 X 10"''. Numbers in parentheses are uncertainty estimates 
for the threshold extrapolation. 
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8 


12 


16 


Eo 


-46.69 


-46.69 


-46.69 


IT-NCSM(l) 


-95.10(2) 


-103.24(2) 


-107.81(2) 


IT-NCSM(2) 


-104.18(15) 


-116.32(15) 


-122.81(50) 


IT-NCSM(2)+MRD 


-104.75(15) 


-117.22(15) 


-123.75(50) 


IT-NCSM(2)+PT(3) 


-104.81(15) 


-117.62(15) 




IT-NCSM(seq) - (a) 


-104.49(10) 


-116.86(25) 


-123.14(70) 


IT-NCSM(seq) - (b) 


-104.43(10) 


-117.12(25) 


-123.45(70) 


full NCSM 


-104.75 







[Hta) and (b). The contribution of the MRD correction grows 
slightly with A^max and reaches about 1 MeV for A^max = 18. 
As seen from Tab.Hl the IT-NCSM(2)h-MRD energy is in ex- 
cellent agreement with the full NCSM. A computationally 
more demanding a posteriori correction based on the explicit 
calculation of the second-order energy contribution on top of 
the IT-NCSM(2) eigenstate as defined by Eq. ([18]) yields very 
similar results. The IT-NCSM(2)+PT(3) energies shown in 
Tab. Uagree very well with both, IT-NCSM(2)-hMRD and full 
NCSM. 

The good agreement with the full NCSM energies is yet an- 
other indication of the efficiency of the importance truncation 
scheme in selecting the relevant configurations and of the reli- 
ability of the threshold extrapolation for recovering the contri- 
bution of excluded configurations. The importance truncated 
space is substantially smaller than full NCSM space as seen in 
Fig. [8jc). For A^max - 8 the importance truncation reduces the 
dimension by two orders of magnitude, for A^max = 12 already 
by four orders of magnitude. This dramatic reduction allows 
us to go to much larger values of A^max than ever possible in 
the full NCSM. 

We can improve the efficiency even further by using the 
sequential IT-NCSM(seq) scheme, which requires only one 
importance update for each value of A^max since it uses a ref- 
erence state constructed from the eigenstate in the A^max - 2 
space. In this way all npn/i-states that are possible in a given 
NmsLxfi^ space are generated in the limit (ATmin, Cmin) — » 0. 

The results of IT-NCSM(seq) calculations for the ground- 
state energy of '^O for = 22 MeV and the sequence of 
A^max values starting from A^max = up to A^max - 22 are pre- 
sented in Fig. |9] We study two different values of Cmin, the 
threshold used in the definition of the reference state, since 
this is the only parameter left after the A-min extrapolation. 
The set of A-min values used for the threshold extrapolation is 
the same as before. 

We observe an excellent agreement with the full NCSM and 
with the IT-NCSM(2) of Fig. [8] The numerical results in 
Tab. Ureveal that the IT-NCSM(seq) energies are slightly but 
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FIG. 9: (color online) Ground-state energy, (a) and (b), and model- 
space dimension (c) as function of A^^ax for ""O obtained within the 
IT-NCSM(seq) scheme using the Vucom interaction for TiQ. = 22 
MeV. Two different values of the parent threshold were used: C„,in = 
5 X 10"* (•) and Cmin = 3 x 10""* (♦). For comparison the results of 
full NCSM calculations with the same Hamiltonian are included (+). 



systematically below the IT-NCSM(2) results. This is due to 
the presence of 5p5h and 6/:'6/j configurations in the model 
space of the IT-NCSM(seq), which are excluded from the IT- 
NCSM(2) space. States beyond the 6p6h level are suppressed 
by the importance truncation, i.e. they do not have importance 
measures above the smallest threshold A-min = 3 x 10"^ used in 
this calculation. The IT-NCSM(seq) calculations for the two 
different reference thresholds Cmax agree within the uncertain- 
ties of the A-min extrapolation, which indicates that the values 
chosen here are sufficiently small to capture all relevant com- 
ponents of the reference state. 

The quality of the IT-NCSM(seq) in comparison to the full 
NCSM is independent of the oscillator frequency fiQ. of the 
underlying basis. As shown in Fig. [TO]both sets of calcu- 
lations are essentially on top of each other The maximum 
deviations are around 300 keV, with the IT-NCSM tending to 
higher energies due to its variational character. 

Based on the results of Fig. |9]we can attempt an extrapo- 
lation A^max °°- Close inspection of the A^max dependence 
reveals a non-exponential behavior for large A^max which af- 
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FIG. 10: (color online) Ground-state energies of obtained for the 
VucoM interaction as function of the oscillator frequency TiQ. for dif- 
ferent Ny„.^^1i£l model spaces. Shown are the results of IT-NCSM(seq) 
calculations (solid symbols) in comparison to full NCSM calcula- 
tions (crosses). 



fects the quality of the extrapolation. This is a property of 
the VucoM interaction used here and is not related to the IT- 
NCSM itself. Similar calculations with other interactions, e.g. 
the chiral N3LO potential after an Similarity Renormalization 
Group evolution used in Ref. [jsl [s^l, do not have this prob- 
lem. If we, nevertheless, use the energies for five consecu- 
tive values of A^max to perform an exponential extrapolation, 
the extrapolated energy has a sizable dependence on the cho- 
sen window in A^max- When using the IT-NCSM(seq) energies 
in the window 14 < A^max < 22 we obtain -133.1 MeV, for 
the range 12 < A^max < 20 we obtain -132.4 MeV, and for 
10 < A^niax < 18 we get -130.8 MeV. In order to arrive at 
a stable extrapolation for the Vucom interaction, one would 
have to go to even larger A^niax or use effective model space 
interactions constructed via a Lee-Suzuki transformation. 



C. Center-of-Mass Contamination 

An important advantage of the NCSM is the possibility 
to exactly separate the intrinsic and the center-of-mass (CM) 
component of the many-body states. Only in this way a non- 
spurious description of the translationally-invariant intrinsic 
state of the nucleus — and all the observables derived from it — 
is guaranteed. As discussed in Sec. IIII Al this property relies 
on the use of a complete Nmaxfi^ model space constructed 
from a harmonic oscillator single-particle basis. Any other 
model-space truncation will destroy the formal separability 
and lead to CM contaminations of the intrinsic states. 

Since the importance truncation reduces the model space 
to a subset of the full Nmaxfi^ space, it might induce a cou- 
pling between intrinsic and CM motion and destroy the exact 
separability. We have to check explicitly that the IT-NCSM 
eigenstates still exhibit the separation between intrinsic and 
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FIG. 11: (color online) Intrinsic ground-state energy of obtained 
in the IT-NCSM(seq) at HCl = 22 MeV using the modified Hamil- 
tonian H/j. Panel (a) shows the intrinsic energies for j8 = (•) and 
/3 = 10 (♦) in comparison to the full NCSM (+). Panel (b) depicts the 
energy difference SEm = £'int(/J = 10) - iJintCS = 0). 

CM motion. 

A well-known tool to probe the presence and extent of the 
coupling is an artificial shift of the excitation spectrum of the 
CM component of the many-body states. Following Gloeck- 
ner and Lawson |57] this can be implemented by adding a 
harmonic-oscillator Hamiltonian with respect to the CM posi- 
tion and the CM momentum Pcm 

1 , mAQ? -, 3 
^™ = 2^""- ^ ^""^ - 2^" ■ ^^'^ 
The modified Hamiltonian 

Hp ^ Hi,, + PH,^ (24) 

is then used instead of the intrinsic Hamiltonian (|22]) at all 
stages of the calculation. 

If intrinsic and center-of-mass motion are properly decou- 
pled, then this shift will not affect the intrinsic state whatso- 
ever. The intrinsic ground-state energy, defined via the expec- 
tation value Eint(J3) - (^F^l //im l^/i) computed with the eigen- 
states 1^^) obtained for H/j, has to be completely independent 
of j6. Any dependence of Ei,t(J3) on /3 signifies an unphysical 
coupling of the intrinsic state to the CM state of the nucleus. 

As an example for this check, we discuss the '^O ground- 
state energy obtained in the IT-NCSM(seq) scheme. The IT- 
NCSM(seq) is set up as described in Sec. lIV Bl Fig.fTTIdepicts 
the intrinsic energies for a sequence of A^n,ax-values obtained 
for y6 = 0, i.e. with the intrinsic Hamiltonian used in all previ- 
ous calculations, and for yS = 10. The intrinsic energies of both 
calculations agree almost perfectly. As shown in Fig. [TTT b) 
the difference is always below 300 keV and consistent with 
within the uncertainty of the threshold extrapolation. Evi- 
dently, the importance truncation does not induce any notice- 
able coupling between intrinsic and CM degrees of freedom 
and thus the eigenstates are free of CM contaminations. 
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The situation is completely different if we start from a 
model space which is not based on the N^^^fiQ. truncation. 
A well known example is the core-plus-valence-space shell 
model, where the model space is spanned by Slater deter- 
minants generated by all possible occupations of a few va- 
lence orbitals. A number of studies show the severity of the 
problem: As discussed in Ref. jssi], e.g., spurious admixtures 
cause the ground-state energy of '''O to be overestimated by 
several MeV. Similar effects are observed when using the im- 
portance truncation idea with a no-core model space defined 
through a truncation of the single-particle basis. These IT-Cl 
calculations, as discussed in Ref. [SD, exhibit sizable CM con- 
taminations of the intrinsic states which also lead to energy 
shifts of several MeV for the '^O ground state. A detailed 
investigation of the CM contaminations in IT-CI and coupled- 
cluster calculations will be presented elsewhere ll59ll . 



V. APPLICATIONS & BENCHMARKS: NON-MAGIC 
NUCLEI 

The IT-NCSM is not limited to doubly-magic or closed- 
shell nuclei. We can apply the same ideas and computational 
techniques, in particular the IT-NCSM(seq) scheme, without 
any changes to non-magic or open-shell nuclei. In this section 
we demonstrate this flexibility and discuss the performance of 
the IT-NCSM scheme for selected non-magic even-even nu- 
clei from the p-shell in comparison to the full NCSM. A sys- 
tematic study of p-shell nuclei with different unitarily trans- 
formed realistic interactions will be presented in a forthcom- 
ing publication. 



A. Carbon-12 

As a first step towards open-shell nuclei we consider the 
ground state of '^C in the IT-NCSM. The computational com- 
plexity of this problem is similar to the "'O ground state, be- 
cause of the incomplete filling of the p-shell. The full NCSM 
is typically limited to A^max = 8, whereas the IT-NCSM can be 
extended to A^max = 22 and beyond. 

Both schemes for constructing the importance truncated 
space, the iterative 1T-NCSM(/) and the sequential IT- 
NCSM(seq) scheme, can be applied without change. For the 
1T-NCSM(;) scheme a natural choice for the initial reference 
state is the ground state obtained from a OfiQ calculation in the 
full NCSM instead of the single Slater determinant that spans 
the OhO. space for a magic nucleus. For the IT-NCSM(seq) 
scheme we start with a full NCSM calculation in a QftQ. or 
2HQ, space in any case, so there is no technical difference be- 
tween closed- and open-shell nuclei. For brevity, we restrict 
ourselves to the IT-NCSM(seq) scheme in this section. As 
in Sec. |IV] we employ a set of calculations with importance 
thresholds in the range /Cmin - 3 x 10"^ to 14 x 10"^ for 
each A^niax- On this basis we perform a constrained thresh- 
old extrapolation as described in Sec. IlllFl making use of the 
second-order perturbative estimate of the energy contribution 
of excluded configurations. 
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FIG. 12: (color online) Ground-state energy, (a) and (b), and model- 
space dimension (c) as function of A'n,„ for '^C obtained within the 
IT-NCSM(seq) scheme using the Vucom interaction for hCl = 24 
MeV. Two different values of the parent threshold were used: Cn,in = 
5 X 10"* (•) and Cmin = 3 x 10"'* (♦). For comparison the results of 
full NCSM calculations with the same Hamiltonian are included (+). 



The evolution of the ground-state energy and the model- 
space dimension with A^max obtained in the IT-NCSM(seq) for 
the Vucom interaction is depicted in Fig. [12] For the refer- 
ence threshold we use two different values, Cmin = 3 x 10""* 
and 5 x 10""*. The sensitivity of the ground-state energy to 
the reference threshold is slightly larger than for the doubly- 
magic because of the absence of a single dominant basis 
state. However, the difference between the two sets of ener- 
gies remains well below 1 MeV. 

As for the general rate of convergence is rather slow 
and of non-exponential character for model spaces beyond 
Mnax ~ 14. To a large extend this can be traced back to the 
high-momentum behavior of the first-generation Vucom inter- 
action. A rough extrapolation based on the five data points 
in the range 14 < A^max ^ 22 leads to an estimated ground- 
state energy of -84.6(1.5) MeV, where the uncertainty is de- 
termined by comparing with extrapolations for other sets of 
five consecutive points. This is almost 8 MeV above the ex- 
perimental ground state energy of -92.16 MeV ll60ll . Keeping 
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FIG. 13: (color online) Ground-state energies of *He (a) and '^He 
(b) as function of A'^ax obtained within the IT-NCSM(seq) scheme 
(•) using the Vucom interaction for TiSi = 24 MeV. For comparison 
the results of full NCSM calculations with the same Hamiltonian are 
included (+). 

in mind that the calculated ground-state energy of '^O was 
at least 5 MeV below the experimental value, this can be in- 
terpreted as evidence for deficiencies in the spin-orbit part of 
the first generation Vucom interactions, which in turn could be 
related to missing three-body interactions. 

Interestingly, a similar pattern has been observed for the 
JISP16 interaction in the full NCSM calculations presented in 
Ref. [6]. Although these NCSM calculations were limited 
to A^max ^ 8 the softness of the JISP16 interaction allows for 
quantitative conclusions already in these small spaces. Based 
on systematic extrapolations the authors conclude that '^C is 
overbound by approximately 2 MeV and '^O is overbound by 
15 to 18 MeV. Hence the difference in the binding energies of 
the two nuclei is of the same order as for the Vucom interaction 
although the JISP16 overbinds '^O significantly. 



B. Helium-6 and Helium-8 

As a second example we consider the neutron-rich Helium 
isotopes ^He and ^He. Whereas for ^He one is able to reach 
large A^^ax with the full NCSM akeady, the few additional 
neutrons in these isotopes significantly reduce the range of the 
full NCSM, typically to A^^ax < 16 for *'He and A^^ax < 12 for 
^He ll6lll . With the importance truncation we can overcome 
this limitation easily. 

The IT-NCSM(seq) results for the ground states of ^He and 
^He obtained with Vucom at HQ. - 24 MeV with Cmm - 
5 X 10 are summarized in Fig. [13] As before, the IT- 
NCSM(seq) energies show an excellent agreement with the 



results of full NCSM calculations where the latter are feasible. 
For larger A^max the threshold extrapolation shows uncertain- 
ties of up to 700 keV for ^He. If necessary, these uncertain- 
ties can be reduced by considering lower /Cmin -values for the 
threshold extrapolation. The general convergence as a func- 
tion of A^niax is rather slow, particularly for ^He. In addition 
to the properties of the Vucom interaction discussed before, 
the structure of these nuclei affects the convergence rate. Ob- 
viously, the description of the neutron halo in an oscillator 
basis requires high-lying single-particle states and thus large 
A^max- Only through the importance truncation these large 
model spaces are accessible. 

Because of the slow convergence and the relatively large 
uncertainties of the threshold extrapolation, an extrapola- 
tion to A^max °° only provides a rough estimate. Using 
the results for the five largest spaces we obtain an extrap- 
olated ground-state energy of -27.4(1.0) MeV for *He and 
of -26.5(1.5) MeV for ^He. A systematic study including 
a variation of the oscillator frequency is needed to provide 
more precise extrapolations. The comparison of these esti- 
mates with the experimental binding energies of -29.27 MeV 
and -31.41 MeV S for ''He and ^He, respectively, con- 
firms our observations regarding the deficiencies of the first- 
generation Vucom interactions. The systematic underbind- 
ing of these open-shell systems could be remedied, e.g., by 
a stronger spin-orbit component of the interaction. Again, the 
NCSM studies with the JISP16 interactions presented in Ref. 
10] show a similar trend, though the absolute deviations are 
smaller. 



VI. CONCLUSIONS & OUTLOOK 

We have introduced an importance truncation scheme with 
all its technical aspects as a new tool to facilitate ab initio 
nuclear structure calculations beyond the domain of conven- 
tional CI approaches. Based on an a priori importance mea- 
sure derived from multiconfigurational perturbation theory we 
identify the important configurations for the description of in- 
dividual target states such that the dimension of the eigenvalue 
problem that needs to be solved is dramatically reduced. The 
effect of excluded configurations can be reliably included by 
combining a perturbative estimate of their energy contribution 
with threshold extrapolation techniques. 

In combination with the Nmaxfi^ space of the NCSM the 
importance truncation provides a powerful tool to asses all 
aspects of nuclear structure in light and medium-heavy nu- 
clei. The importance truncation preserves a crucial property 
of the NCSM: the decoupling of intrinsic and center-of-mass 
degrees of freedom which guarantees that the intrinsic observ- 
ables are free of unphysical center-of-mass contaminations. 
We have discussed two schemes for setting up the importance- 
truncated space, the iterative IT-NCSM(;) and the sequential 
IT-NCSM(seq) scheme. The latter is most efficient since we 
have to construct the importance-truncated space only once 
for each A^n,ax- 

Moreover, the IT-NCSM(seq) scheme is conceptually supe- 
rior, because in the limit of vanishing thresholds C^in and K^i^ 
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the complete Nmax^^ space is obtained without any truncation 
regarding the npnp order at each step of the sequence of A^max 
values. Hence, the full NCSM results are recovered in the 
limit (Cmin, /fmin) ^ at cach Nimx- Based on this property we 
use a numerical a posteriori threshold-extrapolation to obtain 
an approximation to the full NCSM with well-defined error 
bounds. The stability of this extrapolation is greatly enhanced 
by using information on the contribution of excluded config- 
urations from perturbation theory. Further improvements of 
these extrapolation techniques, e.g. along the lines discussed 
in Refs. ll39ll or ll62ll . will be investigated in the future. 

Our series of benchmark calculations confirms the excel- 
lent agreement of the IT-NCSM with the full NCSM in all 
cases where the latter is computationally feasible. The com- 
parison also demonstrates that the IT-NCSM gives access to 
much larger Nmnxh^l spaces and to heavier nuclei than the full 
NCSM. The range of the IT-NCSM in both, N^^^ and A is only 
limited by the computing time and not by memory. Moreover, 
the time-consuming steps of the computation can be easily 
parallelized with minimal communication overhead. 

The present calculations also allow for a detailed as- 
sessment of the first-generation Vucom interactions used. 
Whereas the '*He binding energy is in agreement with experi- 
ment by construction, the ground state of '^O is overboundby 
at least 5 MeV. This level of agreement is still satisfactory and 
is not found with most other realistic two-body interactions, 
be it bare of effective. For the non-magic nuclei discussed 
here the binding energies are systematically underestimated 
with the Vucom interaction, which might hint at deficiencies 
in the spin-orbit part of the interaction. Furthermore, the IT- 
NCSM calculations show that the convergence rate of the first- 
generation Vucom when going to large spaces is rather slow, 
which might result from the high-momentum behavior of the 
interaction. All of these deficiencies will be addressed during 
the construction of the next generation of UCOM-transformed 
interactions and the IT-NCSM provides a indispensable tool 
for assessing these aspects. 

Obviously, the investigation of ground states of closed- and 
open-shell nuclei is only a first step towards a complete ab ini- 



tio description of nuclear structure. The next crucial step is the 
extension of the IT-NCSM to excited states. The importance- 
truncation scheme can be generalized in a straight-forward 
manner for the simultaneous description of a few target states. 
In this way it becomes possible to describe, e.g., ground and a 
few excited states simultaneously and on the same footing. A 
detailed discussion of the methodical details will be presented 
in a subsequent paper, together with a variety of applications. 

Since we automatically obtain a representation of the eigen- 
states in a shell-model basis, all observables of interest can 
be computed directly. Although we discussed only energies 
for the purpose of the present benchmark, we have computed 
a variety of properties, e.g., radii, density distributions, and 
form factors. We have even used the IT-NCSM eigenstates 
as input for the calculation of phase-shifts for low-energy 
nucleon-nucleus scattering reactions in the framework of the 
NCSM/resonating group method (NCSM/RGM) SIH. 

This demonstrates that the IT-NCSM off'ers the same possi- 
bilities for complete ab initio calculations of nuclear structure, 
spectroscopy, and reactions as the full NCSM. At the same 
time, the IT-NCSM extends the range of these ab initio stud- 
ies to heavier nuclei and larger model spaces, which is crucial 
for developing a consistent framework for nuclear structure 
theory throughout the whole nuclear chart. 
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